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A.   Introduction 

The  Loran  system  is  a  radio  aid  to  navigation  which 
utilizes  the  principle  of  hyperbolic  fixing.   The  locus  of 
points  for  which  the  difference  in  arrival  time  of  synchron- 
ized signals  from  a  pair  of  transmitters  is  constant  det- 
ermines a  hyperbolic  line  of  positions  (LOP) .   The 
intersection  of  two  hyperbolic  lines  of  position  from  two 
pairs  of  transmitters  (stations)  determines  a  position 
estimate  (a  fix) .   That  two  lines  of  position  and  hence 
two  pairs  of  stations  are  required  for  a  fix  does  not 
necessarily  mean  that  four  separated  stations  are  required, 
since  one  station  of  one  pair  may  be  colocated  with  one 
station  of  the  other  pair.   In  this  case,  a  Loran    triplet 
is  formed  (Figure  1) .   Triplets  may  be  joined  "end-to-end" 
by  station  colocation  to  form  a  Loran    chain    (Figure  2) . 
Loran  chains  are  common  on  both  the  East  and  West  Coasts 
of  the  North  American  continent. 

"Standard  Loran"  or  "Loran-A"  operates  at  a  fre- 
quency just  below  2  MHz  and  is  now  only  used  in  the 
Pacific  area.   "Loran-C"  operates  at  100-kHz  and  is 
used  both  by  ships  and  aircraft  world  wide.   The  compu- 
tational algorithm  and  programs  described  herein  can  be 
used  for  position  determination  with  Loran-C  triplets. 
Further  information  on  the  history,  development  and 
operation  of  the  Loran  systems  may  be  found  in  References 
1  and  2 . 
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MrM2 


(a)   Colocated  Master  Stations 


M. 


M. 


(b)   Colocated  Slave  Stations 


M, 


(c)   Colocated  Master  and  Slave 


Figure  1.   Loran  Triplets 
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M1'M2 


M3,M4 


S4'S5 


Figure  2.   Loran  Chain  of  Five  Loran  Triplets 


B.   Program  Description 

The  Position  Determination  Program  requires  opera- 
tional data  cards  (described  below)  for  on-location  position 
determination  from  Loran  triplet  time-difference  measure- 
ments.  Two  programs  are  used  to  prepare  an  operational 
data  card.   Operational  data  cards  should  be  prepared  and 
validated  prior  to  on-location  navigational  use.   Then, 
although  three  programs  are  used  in  total,  only  one 
program  is  required  for  position  determination,  the 
remaining  two  programs  are  used  to  prepare  operational 
data  cards . 

The  programs  are : 
Master  Data  Card  Program;   This  program  is  used  to  prepare 
master   data    cards.      A  master  data  card  requires  the  follow- 
ing information  for  a  master  (M)  station/slave  (S)  station 
pair  (a  list  of  current  station  parameters  is  in  the 
appendix) : 

1.  A  M/S  pair  identification  number. 

2.  The  quantity  At  which  is  the  sum  of  the  coding 

delay  plus  the  one  way  base  line  time  in  microseconds. 

3.  The  latitude  and  longitude  of  the  master  station. 

4.  The  latitude  and  longitude  of  the  slave  station. 
Some  preprocessing  of  these  data  is  performed  before  the  master 
data  card  is  generated.   The  data  generated  require  only  one 
side  of  a  TI-59  magnetic  card  for  each  M/S  pair,  thus  a 
second  M/S  pair  may  be  placed  on  side  2  of  a  master  data 

card.   A  master  data  card  should  be  prepared  in  advance  for 
each  M/S  pair  that  might  be  used  within  an  area  of  operation. 


Operational  Data  Card  Program;   This  program  is  used  to 
prepare  an  operational    data    card    for  each  Loran  triplet 
that  might  be  used  within  an  area  of  operation.   Each 
operational  data  card  contains  data  merged  from  two 
master  data  cards  that  contain  the  M/S  pair  information 
for  each  pair  of  a  triplet.   These  merged  data  are  validity 
checked,  colocation  of  master  or  slave  determined  and 
encoded. 

The  inputs  required  to  generate  an  operational  data 
card  for  a  Loran  triplet  using  this  program  are  those 
contained  in  the  two  master  data  cards  that  are  associated 
with  the  triplet.   It  is  possible  to  prepare  and  store 
operational  data  cards  rather  than  master  data  cards.   This 
may  be  desirable  if  there  is  no  scarcity  of  cards  and  storage 
space,  however  the  number  of  possible  Loran  triplets  is 
considerably  larger  than  the  number  of  M/S  pairs. 

Position  Determination  Program:   This  program  is  used  for 
position  determination.   The  required  inputs  for  the  program 
are  the  indicated  time  differences   T   for  the  two  M/S 
pairs  of  the  triplet  and  the  data  contained  in  the  opera- 
tional data  card  associated  with  the  triplet.   The  outputs 
of  the  program  are  the  computed  latitude  and  longitude  of 
the  Loran  fix.   Note:   Every  Loran  determination  gives 
two  solutions:   the  proper  solution  and  the  spurious  solution 


The  spurious  solution  can  almost  always  be  rejected  by 
inspection.   However,  if  the  stations  of  the  Loran  triplet 
are  nearly  aligned  this  may  not  be  the  case,  and  additional 
information  will  be  required  to  determine  the  proper 
solution. 


C.   TI-59  Calculator  Programs 
1.   User  Instructions 


Master  Data  Card  Program  (MDCP) 


Step 


Instructions 


Input 
Data/Units   Keys 


Output 
Data/Units 


ID 

2nd, A' 

ID 

At 

2nd,B' 

At 

*.. 

A 

0 

M 


0.  Turn  the  calculator  on.   If  the  cal- 
culator has  been  on,  turn  it  off, 
then  on  to  initialize  it. 

1.  Read  side  1  of  the  MDCP  program  card. 

2.  Press  CLR  and  then  read  side  2  of  the 
MDCP  program.  CLR 

3.  Input  a  unique  ID  number  for  the  Loran 
pair.   See  Note  a. 

4.  Input  the  coding  delay  At. 

5.  Input  the  master  station  latitude. 
See  Note  b. 

6.  Input  the  master  station  longitude. 
See  Note  b. 

7.  Input  the  slave  station  latitude.  (j)        C  0 

8.  Input  the  slave  station  longitude.  A        D  0 

9.  Run.   When  the  calculations  are  E 
complete,  4  will  appear  in  the 

display;  then  the  display  will  blank. 
After  the  display  has  blanked,  enter 
a  card  in  the  read/write  slot  to 
record  the  master  data  card  for  the 
Loran  pair.   When  the  recording  is 
complete,  4  will  reappear  in  the 
display.   See  Note  c. 

Notes:   a.   Loran  pairs  are  coded  on  the  navigation  maps  using 
designators  such  as  9930X,  9930Y,  9930Z  and  9930W. 
It  is  suggested  that  the  ID's  for  these  pairs  be 
coded  as  9930.1,  9930.2,  9930.3  and  9930.4  respectively, 
However,  any  consistent  scheme  is  acceptable. 


b.  The  format  for  latitude  and  longitude  data  input  is 
of  the  form:   DDD.MMSSFF,  where 

DDD  denotes  degrees 

MMM  denotes  minutes 

SS  denotes  seconds 

FF  denotes  hundredths  of  a  second. 

A  minus  sign  (-)  denotes  South  latitude  and  West 
longitude.   See  the  Appendix  with  regard  to  the 
latter  sign  convention. 

c.  When  the  display  has  blanked  after  4  has  been  displayed, 
the  calculator  is  prepared  to  write  a  card  side  with 
the  information  (data)  contained  in  bank  4.   When  a 
card  is  inserted,  it  will  be  recorded.   Do  not  restrict 
the  cards  motion  once  the  card  has  been  grasped  by  the 
read/write  mechanism.   When  recording  is  complete,  the 
card  will  stop  and  4  will  reappear  in  the  display. 

The  card  can  then  be  removed. 


Operational  Data  Card  Program  (ODCP) 


Step 


Instructions 


Input 
Data/Units 


Output 
Keys   Data/Units 


0.  Turn  the  calculator  on.   If  the  calculator 
has  been  on,  turn  it  off,  then  on  to 
initialize  it. 

1.  Read  side  1  of  the  ODCP  program  card. 

2.  Press  CLR  and  then  read  side  2  of  the  ODCP 
program  card. 

3.  Press  CLR  and  then  read  a  master  data  card 
containing  the  first  pair  of  a  Loran  triplet 

4.  Press  A. 

5.  Press  CLR  and  then  read  a  master  data  card 
containing  the  second  pair  of  a  Loran 
triplet. 


CLR 

CLR 

A 

CLR 


Press  B. 


See  Note  a. 


Run.   When  the  calculations  are  complete,  4 
will  appear  in  the  display;  then  the  display 
will  blank.   After  the  display  has  blanked, 
enter  a  card  in  the  read/write  slot  to 
record  an  operational  data  card.   When  re- 
cording is  complete,  4  will  reappear  in  the 
display.  See  Note  b. 


Notes:   a.  Should   the  display  flash,  then  the  two  master  data 
cards  do  not  correspond  to  a  Loran  triplet.   Both 
the  latitude  and  longitude  of  the  colocated  stations 
must  be  identical  on  both  master  data  cards  in  order 
to  successfully  produce  an  operational  data  card. 

b.  When  the  display  has  blanked  after  4  has  been  displayed, 
the  calculator  is  prepared  to  write  a  card  side  with 
the  information  (data)  contained  in  bank  4;  and,  when 
a  card  is  inserted,  it  will  be  recorded.   Do  not 
restrict  the  cards  motion  once  the  card  has  been 
grasped  by  the  read/write  mechanism.   When  recording 
is  complete,  the  card  will  stop  and  4  will  reappear 
in  the  display.   The  card  can  then  be  removed. 


Position  Determination  Program  (PDP) 


step 


Instructions 


Input 
Data/Units   Keys 


0.  Turn  the  calculator  on.   If  the  cal- 
culator has  been  on,  then  turn  it  off, 
then  on  to  initialize  it. 

1.  Read  side  1  of  the  PDP  program  card. 

2.  Press  CLR  and  then  read  side  2  of  the 
PDP  program  card. 

3.  Press  C  to  initialize,  and  then  read  the 
operational  data  card  for  a  Loran  tri- 
plet that  you  are  receiving.   See  Note 

a  below. 

4.  Input  the  observed  time  delay  from  the 
first  Loran  pair. 

5.  Input  the  observed  time  delay  from  the 
second  Loran  pair. 

6.  Compute  the  1st  solution  latitude. 

7.  Compute  the  1st  solution  longitude. 

8.  Compute  the  2nd  solution  latitude. 

9.  Compute  the  2nd  solution  longitude. 
See  Note  b. 


CLR 


Output 
Data/Units 


D 

latitude 

R/S 

Longitude 

E 

Latitude 

R/S 

Longitude 

10.     To  compute  new  solutions  for  new  time 
delays,  go  to  step  4.   To  compute  new 
solutions  for  a  new  triplet,  go  to 
step  3. 

Note  a.   If  the  data  card  does  not  read  (a  flashing  display  results) , 
press  CLR  and  repeat  the  step.   If  the  card  does  not  read 
on  the  2nd  try,  the  card  should  be  cleaned  (see  the  TI-59 
owners  manual)  or  a  new  data  card  prepared. 

Note  b.   Once  the  proper  solution  (solution  1  or  solution  2)  has 
been  determined  for  a  triplet,  then  only  step  6  (for 
solution  1)  or  step  8  (for  solution  2)  need  be  repeated 
after  inputing  new  time  delays  in  step  4  and  5. 
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2.  Example  Problem  for  Loran-C  Pairs  9930X  and  9930Y 

An  example  use  of  the  Master  Data  Card  Program  using  coding  delays 
and  latitude  and  longitude  values  from  the  Appendix: 


Step 


Instructions 


Input 
Data/Units 


Keys 


Output 
Data/Units 


Turn  the  calculator  on.   If  the  cal- 
culator has  been  on,  turn  it  off,  then 
on.   This  assures  that  it  will  be  in 
the  correct  initial  state . 


Read  side  1  of  the  MDC  program  card. 

Press  CLR,  then  read  side  2  of  the  MDC 
program  card. 


CLR 


Input  the  ID  for  9930X.   Next  press 
2nd,  then  A. 


9930.1 


2nd, A'    9930.1 


Input  the  coding  delay  At  for 
9930X.   Next  press  2nd,  then  B. 


36389.66 


2nd,B'    36389.66 


Input  the  master  station  latitude , 
then  press  A. 


34.034604 


A 


Input  the  master  station  longitude 
and  press  +/-  for  West,  then  press 
B. 


-77.544676 


Input  the  slave  station  latitude, 
then  press  C. 


46.463218 


Input  the  slave  station  longitude 
and  press  +/-  for  West,  then  press 
D. 


-53.102816 


10. 


Press  E.   After  4  has  appeared  in 
the  display,  enter  an  available 
card  side  in  the  read/write  slot  of 
the  calculator.   After  recording  is 
complete,  4  will  reappear.   Then  re- 
move the  card  and  label  the  recorded 
side  "9300X  MASTER." 

Input  the  ID  for  9930Y,  next  press 
2nd  and  then  A. 


9930.2 


2nd, A'   9930.2 


II 


Step 


Instructions 


Input 
Data/Units 


Output 
Keys    Data/Units 


11.      Input  the  coding  delay  At   for 
9330Y,  next  press  2nd  and  then  B. 


52541.31 


2nd,B'   52541.31 


12.      Input  the  master  station  latitude, 
then  press  A. 


34.034604 


13.  Input  the  master  station  longitude 
and  press  +/-  for  West,  then  press 
B. 


-77.544676 


14.      Input  the  slave  station  latitude, 
then  press  C. 


41.151193 


15.  Input  the  slave  station  longitude, 
and  press  +/-  for  West,  then  press 
D. 

16.  Press  E.   After  4  has  appeared  in 
the  display,  enter  an  available 
card  side  (for  example  the  second 
side  of  the  card  used  in  Step  9)  in 
the  read/write  slot  of  the  cal- 
culator.  After  recording  is 
complete,  4  will  reappear  in  the 
display.   Then  remove  the  card  and 
label  the  recorded  side  "9930Y 
MASTER."   The  recorded  data  cards  (or 
card)  will  be  used  in  the  next 
example. 


-69.583909 


D 

E 
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ts    An  example  use  of  the  Operational  Data  Card  Program: 
-~ 

Input  Output 

Step  Instructions  Data/Units    Keys    Data/Units 

0.  Turn  the  calculator  on.   If  the  cal- 
culator has  been  on,  turn  it  off,  then 
on  to  initialize  it. 

1.  Read  side  1  of  the  ODCP  program  card. 

2.  Press  CLR  and  then  read  side  2  of  the 
ODCP  Program  Card.  CLR 

3.  Press  CLR  and  then  read  the  master 
data  card  for  station  9330X.  CLR 

4.  Press  A.  A 

5.  Press  CLR  and  then  read  the  master 
data  card  for  station  9930Y.  CLR 

6.  Press  B.   See  Note  a  of  the  ODCP 
instructions  on  page  9.  B 

7.  Press  E.   When  the  display  has  E 
blanked  after  4  has  appeared,  insert 
an  available  card  side  in  the  read/ 
write  slot.   After  recording  is 
complete,  4  will  reappear  in  the 
display.   Label  the  recorded  card 
side  "9330X/9930Y  OPERATIONAL  DATA. 
Then,  for  identification,  label  the 
A  key  position  on  the  card  "9330X" 
and  the  B  key  position  "9330Y". 
This  card  will  be  used  as  the  data 
card  in  the  Position  Determination 
Program  example.   The  key  labeling 
relates  to  that  program. 


1  5 


An  example  use  of  the  Position  Determination  Program: 
step 


Instructions 


Input  Output 

Data/Units   Keys   Data/Units 


You  are  in  the  North  Atlantic  and  you 
are  receiving  9930X  and  9930Y  and  you 
wish  to  obtain  a  fix. 

0.  Turn  the  calculator  on.   If  the  cal- 
culator has  been  on,  then  turn  it 
off,  then  on  to  initialize  it. 

1.  Read  side  1  of  the  PDP  program  card. 

2.  Press  CLR  and  then  read  side  2  of  the 
PDP  program  card. 

3.  Press  C  to  initialize,  and  then  read 
the  operational  data  for  the  triplet 
9930X/9930Y  from  the  card  side  re- 
corded in  the  ODCP  example. 

4.  The  indicated  time  delay  is  28800  ys 
from  99 3 OX. 

Input  the  indicated  time  delay,  then     28800 
press  A. 

5.  The  indicated  time  delay  is  49400  ys 
from  9930Y. 

Input  the  indicated  time  delay,  then     49400 
press  B. 

6.  Press  D  to  obtain  the  1st  solution 
latitude. 

7.  Press  R/S  to  obtain  the  1st  solution 
longitude. 

The  first  solution  is:  42°44'57"N  Latitude 

41°07'32"W  Longitude 

8.  Press  E  to  obtain  the  2nd  solution 
latitude. 

9.  Press  R/S  to  obtain  the  2nd  solution 
longitude. 

The  2nd  solution  is:    27°00'07"S  Latitude 

102o27'12"E  Longitude 


CLR 


B 


28800 


49400 


D      42.4457 
R/S    -41.0732 


-27.0007 


R/S    102.2712 
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3 .   Program  Listings 


Master  Data  Card  Program 


0  0  0 

76 

LBL 

0  0 1 

42 

STD 

f !  fl  '-' 

■?  ■"■ 

S  T  # 

i     i— 

0  0  3 

0  0 

0  0 

0  0  4 

97 

nsz 

n  r> tr 

0  0 

0  0 

0  0  fe- 

42 

STD 

ll  0  7 

0  0 

0 

0  0  8 

IN';1 

n  n '-? 

cr  O 

F  T  '::' 

'—  :-' 

t            XII 

i  j  1  0 

04 
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Position  Determination  Program 
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Position  Determination  Program  (continued) 
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Position  Determination  Program  (continued) 
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D.   Loran-C  Fixing  Algorithms 

The  development  of  the  Loran  fixing  algorithms  in  this 
report  is  presented  in  more  detail  in  a  companion  report 
[Ref.  3]  and  will  not  be  repeated  here. 

The  basic  Loran-C  equation  [Ref.  4]  can  be  written  as 

T  =  [Tfi  +  p(TQ)]     [TM  +  P(TJ]  +  [Tn  +  p(TM)]  +  6    (1) 

where 

T   is  the  "indicated  time  difference"  in  microseconds, 
TM,Tq   is  the  distance,  in  microseconds,  from  the  master 
and  the  slave  to  the  receiver,  respectively, 
T_.   is  the  distance,  in  microseconds,  between  the 

master  and  the  slave, 
6    is  the  assigned  coding  delay,  in  microseconds,  and 
p(T)   is  the  secondary  phase  correction,  in  microseconds, 
for  an  all  sea  water  path  of  length   T. 


The  quantity 


At  =  [TB  +  P(TB) ]  +  6 


is  a  constant  for  each  master/slave  pair.   The  following 
World  Geodetic  System  1972  (WGS  72)  values  have  been  adopted 
for  Loran-C  navigation  [Ref.  4]: 
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v„  =  299792458  meters/second  is  the  velocity  of  light 
in  free  space, 

ii   =  1.000  538  is  the  index  of  refraction  of  the  surface 
of  th'i  earth  for  standard  atmosphere  and  100kHz 

t'luct-oraagentic  waves, 

a   -  6378135.00  meters  is  the  equatorial  radius  of  the 
e 

earth,  and 

f   =  1/298.26  is  the  flattening  factor  (1  -  b/a  ,  where 
b   is  the  polar  radius)  of  the  earth. 

Accurate  formulas  for  computing  the  secondary  phase 
correction   p(T)   are  contained  in  Reference  4,  but  for  use 
with  the  handheld  calculator  the  following  linear  approximation 
[Ref.  3]  will  be  used: 

p(T)  =  ax  +  a2T  , 

where 

a,  =  -0.321, 
and  L 

a2  =   0.000635. 

Using  this  approximation,  it  is  possible  to  solve  Equation  1 

for  the  quantity   T0  -  TK4.   We  find 

S     M 

Ts  -  TM  =  (T  -  At)/(1  +  a2)  .  (2) 

On  the  surface  of  a  sphere  a  hyperbolic  line  of  position 
can  be  represented  by  the  equation  [Ref.  3,  page  175] 
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tan  r  =  -, cos  2a  -  cos  2c 

sin  2c  cos  w  +  c  sin  2a 


where  the  origin  of  the  coordinate  system  is  at  the  prime  focus 
of  the  spherical  hyperbola,  2c   is  the  spherical  arc  joining 
the  foci,  2a   is  a  constant  for  any  one  LOP,  and   r   and   lo 
are  the  spherical  coordinates  of  a  point  on  the  LOP.   If  the 
base  line  of  the  coordinate  system  is  the  arc  joining  the  foci 
then   w   is  the  spherical  polar  angle  from  the  base  line  to  a 
point   P   on  the  LOP  and   r   is  the  spherical  polar  distance 
(or  arc)  from  the  prime  focus  to   P.   Using  the  Loran  system 
we  take    f,  =  +1   if  the  prime  focus  is  at  a  master  station 
and  we  take  C  =  -1   if  the  prime  focus  is  at  a  slave  station. 

If  we  take   v  =  vn/n   to  be  the  velocity  of  100kHz  electro- 
magnetic radiation  of  the  earth's  surface  then 


or,  using  Eq .  (2)  , 


where 


2a  =  v(Ts  -  TM)/ae  , 


2a  =  (T  -  At) /a   ,  (4) 


a  (1  +  a  ) 

a   =  _® — £_  _  21295.87  ys  . 

P      vQ/n 

The  baseline  between  master  and  slave  can  be  obtained  from 

2c  =  v  TB/ae  .  (5) 

Here   2c   is  computed  by  program  card  1  (preparation  of  master 
data  cards)  using  the  algorithm  in  Section  E. 
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Consider  a  Loran-C  triplet  with  master  stations  colocated. 
Let  f.         and   '    denote  the  azimuth  angles  of  slave  1  (S..) 
and  slave  2  (S  ),  respectively,  measured  from  North  toward  the 
East  from  the  master  stations  (M)   (see  Fig.  3) .   Further,  let 
u   and   r   denote  the  azimuth  and  spherical  polar  arc  (distance) 
of  the  receiver  (R)  from   M.   For  this  geometry,  Eq .  (3)  can 
be  written  as 

B. 

tan  r.  =  ■= -, = — v — ; — = —  (6) 

l    C .  cos (a-£ . )  +  A . 


where 


A .  =  c,  .    sin  2a . 
li        l 


B.  =  cos  2a.  -  cos  2c. 
ill 


and 


C .  =  sin  2c  . 

l  l 


for  the  i —  Loran  pair,  i  =  1,2.   Since   r  =  r,  =  r~,  tan  r. 

1     2        i 

can  be  eliminated  in  Eq.  (6) .   The  resulting  equation  can  be 
rewritten  as 

C  cos  u  +  S  sin  a  =  K  ,  (7) 


where 


C  =  B,C2  cos  £_  ~  B2Cl  cos  ^l' 
S  -  B1C2  sin  C2  "  B2ci  sin  ^i' 
and      k  =  B2Al  "  BiA2  * 
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N 


M 


Figure  3.   Geometry  of  a  Loran  Triplet  and  a  Receiver 
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J  f  we  define   ()  >  0   and   ,   by  the  equations 


and 


then 


and 


[>    cos  v  -  C, 
[>    sin  ,  =  S, 


(8) 


m  =  ^2  +  s2 , 


y  =  qatn(S,C) 


Here  the  function   qatn(y,x)   is  the  arctangent  of   y/x 
adjusted  for  the  proper  quadrant  according  to  the  signs  of   x 
and   y.   A  compact  form  of  this  function  is 


qatn(y,x)  =  tan 


-1 


x  +  10  9t (x  -  0?) 


+  nt(x  <  0?) 


where 


and 


t(z)  =  1   when   z   is  true 


t ( z)  =  0   when   z   is  false 


When  convenient  we  will  use  the  notation   qatn(y/x)   inter- 
changeably with  qatn(y,x).   The   qatn   function  is  equivalent 
to  the  polar  angle  obtained  using  the  rectangular  to  polar 
conversion  function  on  the  TI-59. 

Now  substitute  Eq .  (8)  into  Eq .  (7)  and  solve  for 


a  =  y  +  cos   (</p) 


9) 
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to  obtain  the  azimuth  angle   ut   of  the  two  points  of  intersection 
of  the  LOP's.   Finally  we  obtain  a  value  for   r   by  substituting 
each   u  into  Eg.  (5) .   We  find  that 


r  =  qatn 


B  . 

1 


C.  cos  (u-f,  .  )  +  A. 
l         l      l 


for   i  =  1  or  2 


The  distance  and  azimuth  from  M  or  the  triplet  vertex  can  be 
converted  into  the  latitude  and  longitude  of  the  two  possible 
positions  of   R. 

The  fixing  algorithm  then  uses  a      and   r   in  the  direct 
solution  algorithm  of  spheroidal  geodesy  (Section  F) . 
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E •   The  Reverse  (Inverse)  Solution  Algorithm 

This  i'i   /'./.-.   solution  algorithm  is  a  modification  of  the 
first  order  in  flattening  (f)  algorithm  given  by  Thomas  [Ref.  5, 
pp.  8-10).   Thomas'  notation  has  been  followed  as  closely  as 
possible  for  ease  of  comparison  of  the  algorithms.   The   qatn 
function  is  defined  in  Section  D.   West  longitudes  (A)  and 
South  latitudes  (()  are  negative.   We  are  given  the  points 
P  ((jj  ,  A ,),  P2((Jj  ,A  )   on  the  spheroid  and  are  to  find  the 
distance   S   between  the  points  and  the  forward  and  back 
azimuths,  u,~   and   u„,  .   Given  quantities  are   4>,,  A  ,  $ 
and   A  .   No  assumptions  about  the  relative  location  of   P, 
and   P2   are  required.   The  modified  reverse    solution  algorithm 
is : 

0   -  tan"  f(l-f)  tai.  <J> .  ]  ,    i  =  1,  2, 

e  -  (o  +e  )/2,  ao  -  (e  -e1)/2,  aa  =  a  -  a1( 

m     i    <s        m     2   1  21 

2  2  2  2 

AA   -  AA/2,   H  =  cos  A0   -  sin   9   =  cos   0   -  sin   A0   =  cos  0   cos  0  , 
m  m         m         m  m        1       2 

2  2  2  ? 

L  =  sin  A0   +  H  sin   AA   =  sin  (d/2),  1  -  L  =  cos  (d/2), 

-1  2        2 

d  =  cos   (1  -  2L),  U  =  2  sin   0   cos   A0  /(l  -  L)  , 

m        m 

2        2 
V  =  2  sin   A0   cos   0  /L,  X  =  u  +  V,  Y  =  U  -  V. 
mm  ' 

T  =  d/sin  d,   6-jd  =  f(TX  -  Y)/4,  S  =  a  (T-6^)  sin  d, 
F  =  2[Y  -  (1  -  2L)  (4  -  X)  ],    G  -  f  T/2 , 
Q  =  -(FG  tan  AA)/4,   AA '  =  (AA  +  Q) /2 
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t      =   qatri(-sin   At)      cos    AA  '  ,    cos    t)      sin   AA  '  )  , 
1  m  m  m  m 


t      =   qdtn(cos   At)      cos    AA  '  ,    sin   t)      sin   AA  '  )     , 
2  in  m  m  m 


°12    ;      S     +    V       «21    =    S    "    S 


This  reverse  solution  algorithm  is  used  by  program  card  1 

(preparation  of  master  data  cards)  to  compute  the  baseline 

distance   2c   and  the  azimuths   £wf,   and   r     between  the 

MS         SM 

master  and  slave  stations  of  a  Loran  pair. 

Details  of  the  modifications  made  to  Thomas'  algorithm 
are  contained  in  Reference  3. 
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F.   The  Direct  Solution  Algorithm 

This  dii'Ci't     solution  algorithm  is  a  modification  of 
the  first  order  in  flattening  (f)  algorithm  given  by  Thomas 
[Ref.  5,  pp.  7-8].   Thomas'  notation  has  been  followed  as 
closely  as  possible  for  ease  of  comparison  of  the  algorithms. 
The   qatn   function  is  defined  in  Section  D.   West  longitudes 
and  South  latitudes  are  negative.   We  are  given  the  point 
P,(({,,,^1)   on  the  spherioid,  where   <J>  ,  *    are  the  geodetic 
latitude  and  longtiude  (geographic  coordinates);  the  forward 
azimuth   a19   and  the  distance   S   to  a  second  point   P_(<j)_,A_)  ; 
and  from  these  we  are  to  find  the  geographic  coordinates 
^0,     \         and  the  back  azimuth   cx9,  .   The  given  quantities  are 
4>  ,  A  ,  ex..  9   and   S.   No  assumptions  about  the  relative  loca- 

X  X  X  £ 

tion    of      P,       and      P~      are    required.       The   modified   direct 
solution    algorithm    is: 

6     -   tan       [(1-f)    tan   41,  ]  .  M  =   cos    0      sin   a 

2 
N    =    cos    0      cos   u  c      =    fM,       c      =    f(l    -   M   )/4, 

D  =    1    -    2c      -   cM,         P   =    c   /D,         0      =   qatn(N,    sin   0    ) 

d   =    S/(a   D) ,         u   =    2(0    -d),         W  =   1    -   2P   cos    u, 

e  1 

V   =    cos(u    +    d)  ,         Y    =    2PVW    sin   d,         Aa    =    d-Y , 
u        =   qatn[-M,    -    (N   cos  Aa    -    sin  0      sin  Au ) ]     , 

K   -    (1-f)  [M2    +    (N   cos  Ao    -    sin  0       sin  Ao) 2] 1/2    , 

4>      =   tan      [  (sin  0      cos  Aa    +   n   sin  Aa  )  /K]    , 

An    -  qatn(sin  Aa    sin  a  cos  6      cos  Aa    -   sin  6      sin  Aa    cos     a      \ 

U   =   CjAo,        AA   -   An   -    h,  A2   -  \y  +  AA. 
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This  direct  solution  algorithm  is  used  by  program 
card  3  (improved  fix  program)  to  compute  the  latitude  and 
longitude  of  the  receiver  using  the  azimuth  and  range  of  the 
receiver  from  the  Loran  triplet  vertex. 

Details  of  the  modifications  made  to  Thomas*  algorithm 
are  contained  in  Reference  3. 
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G.   Some  Results 

The  following  procedure  was  used  to  test  the  fixing 
algorithm:   A  more  accurate  version  of  the  reverse  solution 
algorithm  (see  Reference  3)  was  used  to  compute  time  delays 
at  selected  positions  from  Loran-C  pairs.   These  computed 
time  delays  do  not  depend  upon  a  fixing  algorithm.   The 
fixing  algorithm  was  tested  by  determining  the  distance 
betv/een  the  fix  and  the  positions  used  to  compute  the  time 
delays.   Typical  results  are  presented  in  Tables  I  through 
IV.   As  can  be  seen,  all  differences  (A)  are  within  1  n.mi. 
The  greatest  difference  (A  =  0.92  n.mi.  in  Table  IV)  was 
obtained  at  a  position  located  2270  n.mi.  from  the  triplet 
vertex. 
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Table  I.   9940X/9940Y 


Position 

Calculated 

Time  Delay 

TI-59 

Fix 

A 
n.mi 

Lat 

Long 

9940X 

9940Y 

Lat (N) 

Long (W) 

24°N 

122°W 

27726.19 

40912.76 

23°59,55" 

122°00'01" 

0.08 

26 

122 

27715.97 

40998.39 

25°59'57" 

122°00*01" 

0.05 

28 

122 

27702.41 

41117.84 

28°00,00" 

122°00'00" 

0.00 

30 

122 

27683.53 

41291.85 

29°59'59" 

122°00'00" 

0.02 

32 

122 

27655.47 

41555.46 

32°00'00" 

122°00'00" 

0.00 

34 

122 

27609.63 

41959.57 

34"00'00" 

122°00'00" 

0.00 

36 

122 

27523.56 

42544.11 

36°00'00" 

121059'59" 

0.01 

38 

122 

27334.61 

43248.22 

38°00'00" 

122059'59" 

0.01 

Table  II.   9940Y/9940W 


Position 

Calculated 

Time  Delay 

TI-59 

Fix 

A 

Lat 

Long 

9940Y 

9940W 

Lat(N) 

Long (W) 

n.mi 

37°N 

122°W 

42892.86 

16257.23 

36°59'59" 

122°01'01" 

0.02 

37 

125 

43056.68 

15765.13 

37°00'00" 

125"00,00" 

0.00 

37 

128 

43137.78 

15327.12 

37°00'00" 

128°00'00" 

0.00 

37 

131 

43191.10 

14970.77 

37°00'00" 

131°00'00" 

0.00 

37 

134 

43232.38 

14683.74 

37°00,00" 

134°00'00" 

0.00 

37 

137 

43267.42 

14449.40 

37°00'00" 

137°00'00" 

0.00 

37 

140 

43298.80 

14254.02 

37°00'00" 

140°00'00M 

0.01 

37 

143 

43327.85 

14087.43 

37°00'01" 

142°59,59" 

0.02 
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Table    III.       7930Z/9930X* 


Position 

Calculated 

Time  Delay 

TI-59 

Fix 

A 
n  .mi 

Lat 

Long 

7930Z 

9930X 

Lat (N) 

Long (W) 

60°N 

30°W 

52437.86 

28451.72 

60o00'03" 

29°59'32" 

0.24 

58 

35 

51960.93 

28391.50 

58°00,00" 

34059'48" 

0.09 

56 

40 

50992.37 

28359.15 

55°59'59" 

39059.54.. 

0.06 

54 

45 

49292.46 

28370.85 

53°59,59" 

44°59,57" 

0.03 

52 

50 

47165.60 

28490.64 

52°00'00" 

49059'57" 

0.01 

50 

55 

45236.59 

29070.48 

50°00,00" 

55o00'00" 

0.00 

48 

60 

44505.60 

30991.94 

48°00,00" 

60°00'00" 

0.00 

46 

65 

44475.70 

33697.14 

46°00'00" 

65o00'00M 

0.00 

44 

70 

44588.91 

36567.42 

43o59.59.. 

69°59,59" 

0.02 

Table 

IV.   9930V 

7/9930X* 

Position 

Calculated 

Time  Delay 

TI-59   ] 

?ix 

A 
n.mi 

Lat 

Long 

9930W 

9930X 

Lat (N) 

Long  (W) 

9°N 

47°W 

13058.04 

36466.46 

8°59'59" 

46059'22" 

0.92 

12 

52 

12984.71 

37288.35 

11°59,34" 

51°59'37" 

0.57 

15 

57 

12898.73 

38267.58 

I4059. 45" 

56°59'47" 

0.28 

18 

62 

12793.91 

39431.32 

17°59'52" 

61°59'54" 

0.16 

21 

67 

12656.52 

40794.36 

20°59,56H 

66°59'57" 

0.08 

24 

72 

12451.30 

42330.55 

23°59'59" 

71059'59" 

0.02 

27 

77 

12097.12 

43876.62 

27o00'01" 

77°00'00" 

0.02 

30 

82 

12973.95 

44768.53 

30°00'01" 

82°00'06" 

0.09 

*The    9930    system  was   disestablished    30   September   1979.      It    is    included  here 
for   illustration   only. 
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APPENDIX.   Loran-C  Station  Parameters 

The  following  list  contains  the  pertinent  parameters 
for  each  Loran-C  station  pair.   This  list  was  compiled  from 
data  in  Reference  4.   Each  row  contains  the  following  inform- 
ation: 

1.  The  Loran-C  station  pair  designator 

2.  At,  the  sum  of  the  coding  delay  plus  one  way  base- 
line time,  including  the  secondary  phase  correction 
for  an  all  seawater  path,  in  micrseconds . 

3.  The  master  station  latitude. 

4.  The  master  station  longitude. 

5.  The  slave  station  latitude. 

6.  The  slave  station  longitude. 

In  this  report,  the  convention  that  a  minus  (-)  sign  denotes 
South  latitudes  and  West  longitudes  has  been  followed 
throughout  Some  users  may  wish  to  reverse  the  longitude 
convention  so  that  (-)  longitudes  are  East  rather  than  West. 
This  may  be  accomplished  by  consistently  changing  the  sign 
of  the  longitude  everywhere    that  it  is  required  for  data 
input;  no  program  modifications  are  required.   If  the  sign 
convention  for  longitude  is  reversed,  care  should  be  taken 
that  the  signs  of  all  longitudes  in  the  list  are  reversed. 
In  columns  3  through  6  the  latitudes  and  longitudes  appear 
to  be  in  decimal  form,  but  the  actual  format  is  DDD.MMSSFF 
where 
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DDD  designates  degrees, 

MM  designates  minutes, 

SS  designates  seconds,  and 

FF  designates  hundredths  of  seconds 
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4990  ft 

. 15972 

4990V 

•3435  3 

593uft 

.13131 

■5930V. 

-, ,-.  ~>  cr  cr 

5990ft 

■  13343 

5990V 

•23937 

5990Z 

•  42266. 

793  0W 

1 5063 

7930ft 

27303 

7930Z 

43212 

7960  ft 

1  3&U4 

7960V. 

29651 

7978W. 

3U065 

7970ft. 

15048 

7970V 

43944 

7970Z. 

63216 

7980U 

12303 

793uft 

27443 

~?."i~ii"hI  1 

4  520  I 

79S0Z. 

61542 

7990ft. 

12755 

799U  V 

.'  2  2  f  -• 

7990Z. 

50999 

8970N 

14  355 

3970ft. 

31162 

39  79  V 

4775  :' 

9930U 

13635 

9930'! 

3633"? 

9930V. 

52541 

993 0Z. 

68560 

9940W. 

1  3736 

9940ft. 

23  3i94 

9940V. 

41967 

9960 W. 

13737 

3360ft 

26969 

9960V 

42221 

9960Z. 

57162 

9970 W. 

1528  3 

9970ft. 

36635 

9970V. 

5946  3 

9970Z. 

30746 

9990!-! 

14875 

9990  V 

32069 

9990 Z. 

46533 

.  13 

[6. 

.  88 

46. 

.  93: 

46. 

■  0 1 

51. 

.  3t 

51. 

.  63 

51. 

.  02. 

59. 

■  1*  l' 

59. 

.  20. 

59. 

.45 

6  3 . 

.  14 

6  -' . 

.64 

6  2 . 

.10. 

62 . 

•  J  j 

C'  2  a 

.  30 

6  2 . 

.54. 

30. 

.  38 

30. 

O  O 

30. 

"7  7' 

3  O . 

.  37. 

38. 

.  30 

~T-  C1 

.  69 

38. 

.  11 

'-\ '  -' 

.  06 

3  9 . 

.74 

89 

.5] 

34. 

.  66 

34. 

■"-'  1 

34. 

"7  "' 

34. 

.  90 

3  9 , 

.  50 . 

3  3 . 

.  3u 

33 . 

.  20. 

42. 

9  ~": 

42. 

.65 

42. 

.  06 

42. 

.  94 

24. 

.  1 2 

24. 

.  1  3 . 

24. 

.  73 

24. 

~i~i 

cr  — ■ 

-  03 

cr— • 

.10, 

cr-:1 

444  395 
444395 

-\6272iJ 
•■132720 

cr  -?  cr  .-1  — ■ .-. 
_'  i"  J  O  1  C' 

cr  — >  cr  .-1  — • .-. 
■J  1  JO  1  O 

575378 

531727 

591727 

591727 

194231 

194231 

175363 

175968 

175968 

175963 

593874 

533874 

53  3374 

59  3874 

522061 

522^61 

522061 

510754 

510754 

5 1 0754 

034604 

0  34604 

0  346314 

034  604 

330662 

330662 

330662 

425060 

425060 

425060 

425060 

48041 

48041 

48041 

48341 

090333 

030983 

090^88 


-169. 
-163. 
-U67. 
■067. 
-122. 
-122. 
-122. 
-045. 
-045. 
-045. 
-142. 
-142. 
-007. 
-007. 

"007. 

-007. 

-085. 
-085. 
-085. 
■085. 

016. 

0 1  6 . 

0 1  6 . 
-387. 
-037. 
-087. 
-077. 
-077. 
-377. 
-077. 
-113. 
■US. 

113. 
■^76. 
-076. 
-076. 
■076. 

141. 

141. 

141. 

141. 
-170. 
■170. 

170. 


303120 
303120 

cr  cr  ■-.  —>  -?  i 
J  J  -'  1  1  1 
crcr  -,-?—>  -1 
■J  J  _■  1  1  1 

220224 

220224 

220224 

102747 

102747 

102747 

433190 

483190 

042671 

042671 

042671 

342671 

100930 

1 00930 

100930 

100930 

430596 

4  30536 

430596 

291214 

291214 

291214 

544676 

544676 

544676 

544676 

495637 

495637 

4356  37 

493386 

493386 

493386 

433386 

19290 

19290 

19290 

18290 

145331 

145981 

145931 


20 


41 
46 


47 
50 
64 

62 

46 


54 
66' 
64 
70 
30 
26 

34 

-■cr 
2 '■— 1 

40 
42 
30 
42 
46 

46 
41 
39 
47 


46 
41 
34 
33 
24 
42 
26 
09 
52 

,j-cr 


144916 

234177 

151133 

4632 1 3 

262085 

034799 

362372 

542653 

175968 

463213 

262021 

262035 

432380 

380615 

542653 

545261 

433302 

315501 

015343 

034604 

312088 

532035 

033643 

593374 

425060 

364334 

015843 

463218 

151193 

510754 

034739 

465699 

191313 

432720 

151193 

034604 

510754 

17077 

443700 

362499 

324566 

494505 

144012 

262021 


-178. 
-069. 
-053. 
-131. 
-113. 
-127. 
-023. 

-007. 

-053. 
-152. 
-131. 

+  008. 

+014. 

-023. 
-008. 

-030. 

-037. 

-080. 

-077. 
012. 
027. 

003. 

-085. 
-076. 
-094. 

~080. 

-05  3. 
-063. 
-037. 
-113. 
-122. 
-114. 
-067. 
-069. 
-077. 
-087. 

153. 

143. 

126. 

138. 
+  173. 
-166. 
-152. 


530970 
173020 
533309 
102316 
151965 
443953 
212335 

crcr-.  1  -?cr 

042671 
102816 
221122 
151965 
173633 
274700 
552175 
435363 
434360 

500003 

065352 

544676 

312396 

520152 

12159U 

100330 

433386 

331347 

065352 

102816 

563309 

231214 

443353 

234453 

431743 

553771 

583909 

544676 

291214 

585 1 5 

430906 

085621 

095523 

105231 

531447 

221122 
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Coverage  of  Loran-C  Systems 


Station 

4990 

5930 

5990 

7930 

7960 

7970 

7980 

7990 

8970 

9930 

9940 

9960 

9970 

9990 

Location 
Central  Pacific 
East  Coast,  Canada 
West  Coast,  Canada 
North  Atlantic 
Gulf  of  Alaska 
Norwegian  Sea 
Southeast  U.S.A. 
Mediterranean  Sea 
Great  Lakes 
East  Coast,  U.S.A. 
West  Coast,  U.S.A. 
Northeast  U.S.A. 
Northwest  Pacific 
North  Pacific 
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Position   determina- 
tion with   Loran-C   trip- 
lets  and    the  Texas    in- 
struments  TI-59  pro- 
6    ocr9^™16516   calculator. 
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56O  Shudde 

$$6t  Position  determina- 

tion with  Loran-C  trip- 
lets and  the  Texas  in- 
struments TI-59  pro- 
grammable calculator. 
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Position  determination  with  Loran-C  trip 
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